Spin nematic states in spin-1 antiferromagnets with easy-axis anisotropy 
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It is well known that spin nematic phases can appear in either frustrated magnets or in those 
described by Hamiltonians with large exotic non-Heisenberg terms like biquadratic exchange. We 
show in the present study that non-frustrated spin-l ID, 2D and 3D antiferromagnets with single-ion 
easy-axis anisotropy can show nematic phases in strong magnetic field. For ID case we support our 
analytical results by numerical ones. 
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PACS numbers: 75.10.Jm, 75.10.Kt, 75.10.Pq 

1. Introduction. Frustrated spin systems have of- 
fered in recent years a wealth of opportunities for the 
study of a broad range of novel types of states and 
phase transitions. Spin nematic phases form a class 
of objects in this area which has received much at- 
tention. Spin nematic states are spin-liquid-like states 
which show a multiple-spin ordering without the conven- 
tional long-range magnetic order. The two-spin order- 
ing can be generally described by the tensor'^' Q"; = 

(SfSJ^) - SafsiSjSi)/^. The symmetric part of Q"f 
describes a quadrupolar order which has been exten- 
sively studied both theoretically and experimentally in 
frustrated systems with ferromagnetic (FM) and anti- 
ferromagnetic (AF) nearest-neighbor and next-nearest- 
neighbor couplings, respectively, in strong magnetic field 
h (see, e.g., Ref.'^l and references therein) and in mag- 
nets with large non-Heisenberg spin couplings such as 
biquadratic exchange (8182)^.'^ It has been also shown 
recently that quantum fluctuations accompanied by a siz- 
able single-ion easy-axis anisotropy can also stabilize a 
nematic phase in the kagome spin-1 antifcrromagnet at 

It is well established that the attraction between 
magnons caused by frustration is the origin of quadrupo- 
lar and multipolar phases in quantum magnets.'^ In par- 
ticular, the bottom of the one-magnon band lies above 
the lowest multi-magnon bound state at ft, = /is, where 
hg is the saturation field, as a result of this attraction 
in magnets with FM and AF couplings between nearest- 
and next-nearest neighbors, respectively. Then, transi- 
tions to nematic phases &i h < hg in such systems are 
characterized by a softening of the multi-magnon bound- 
state spectrum rather than the one-magnon spectrum. 

We show in the present paper that the niagnon attrac- 
tion arises also in spin-1 ID, 2D and 3D non-frustrated 
AFs with easy- axis single-ion anisotropy that leads to 
stabilization of nematic phases in strong magnetic field. 
We support our analytical results by numerical ones in 
the particular case of AF chain. 

2. Model and technique. We discuss axially sym- 
metric spin-1 systems described by the Hamiltonian 



H = Y. J.,8,8, +DY,{Stf - hJ2s^ 



(1) 



p + q/2 k + q/2 

r(q,p,k) = \f = 

-p + q/2 -k + q/2 



z + q/2 





-z+q/2 



FIG. 1. Bethe-Salpeter equation for the four-particle vertex. 
Black points stand for bare vertices given by Eq. Jsl. 



where («, j) denote spin pairs coupled with exchange con- 
stants JijO, and D < is the value of the single-ion easy- 
axis anisotropy. 

We examine in the present paper the possibility of the 
nematic phase formation below the saturation field hs by 
considering the transition from the fully polarized state 
that is discussed using the Holstein-Primakoff transfor- 
mation 






(2) 



Expanding the square root in Eq. ([2|, putting all opera- 
tors a\ to the left of all a^ using commutation relations, 
discarding terms containing more than three operators 
a\ and Oi (this is reasonable because one can neglect in- 
teraction of more than two particles in a dilute gas of 
magnons which arises at ft. « ft^), Inland substituting the 
resulting expressions for S~ and S^ into Hamiltonian (fl]) 
one obtains 



•H2 = XI ^^■^P - '^Jo - ^(^-^ - 1) + h] ala 



P' 
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H4 = 



N '^ . 
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-Jl+3-TS{Ji + J3 



a{a2a_3a_4. 
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where N is the number of spins, J- = 1 — \J\ — l/^S, 
the momentum conservation law ^^ p^ = is implied in 
Eq. (pi), we omit some indexes p, and we set Ho = in 
the subsequent discussion. 

At h > hg , the one-magnon spectrum and the magnon 
Green's function are given solely by 'H2, they are exact 



and have the form 



Cp = SJp - SJo - D{2S - 1) 



G(w,p) 



uj — Cp + iS 



(6) 
(7) 



Two-magnon bound states are examined via analysis of 
the pole structure of the two-particle vertex function 
r(a;, q, p, k) (see Fig. [I]), which can be found analytically 
at /i > hg. It can be shown that T{uj, q, p, k) is given by 
a series of ladder diagrams in the fully polarized phase 
which lead to the Bethe-Salpeter equation for the vertex 
shown in Fig. Ill One can represent Jp in the following 
way on any Bravais lattice: 



diUpdLOkGip + q/2)G{-p + q/2)G{k + q/2)G{-k + q/2) 

= 1/{UJ - ep+q/2 - e-p+q/2)(^ - ek+q/2 - e-k+q/2)- 

Consequently, r(a;, q, p, k) has the form 

/,(p)/;(k) 



r(a;-e2(q),q,p,k) 



cj - e2(q) +i& 



(13) 



and the wave function (not normalized) is related to 
/g(p) as follows: 



*,(?) 



hip) 



(14) 



Jp — 2 J" cos Pa, 



(8) 



where a enumerates exchange constants and summation 
over repeated Greek indexes is implied here and below. 
Using Eq. ([8]) we can represent the equation shown in 
Fig. [ijas follows: 

r(c^,q,p,k)=rO(c^,q,p,k) 

D + J°' cos Pa COS Za — 2FSJ°' COS %^ [cos Pa + COS Za) 



e2(q) - ep+q/2 - e-p+q/2 

Certainly, it can be obtained also from Eq. ^ that the 
vertex r(a; ^ e2(q),q, p,k) near the pole can be repre- 
sented in the form (13). It is also evident from Eq. ^ 
that /g(p) has the form 



/?(?) = Aq + 3"^ cos Pa 



(15) 



Substituting Eqs. ( 13 1 and ( 15 ) into Eq. ^ and neglect- 
ing r°(a;, q, p, k) we obtain the following set of equations 



I?(g,z) 



for Aq and Bq-. 



d'^z 

xr(w,q,z,k)^^, 



(9) 



Miq,D,h,iv) 



B 



V{q, z) = 2SJ^ cos ^ cos zp - SJo - D{2S - \) + h - -, 

(10) 

where d is a lattice dimensionality. Let us consider the 
following two-particle Green's function: 



9 
3a=2 



V 



(16) 
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where M(g, D, h, cj) is a square matrix with the following 
elements: 



Moo = l + Dg- 2TSJ" cos 



q 



G"(c.,q,p,k)= r 

Jo 



dt 



-i{H-L0)tJ 



(11) 



Q 
Moa — DVa — 2FSJ^ COS —1 13a, 



J 



^q/2+k"q/2-k 



|0), 



MaO — JaVa — '^Q^'SJa COS 



qa 



(17) 



where 10) denotes the fully polarized state. 
Eq. (11) coincides with the conventional def- 
inition of the two-particle Green's functiorF' 

(r(aq/2+p(ti)aq/2_p(i2)aq/2+k(^3)aq/2_k(*4))> at 

ti — t2 — t and t'i — t^ — 0, where ap(i) denote 
operators in the Heisenberg representation. If there is a 
two-magnon bound state |2;q) with momentum q and 
energy e2(q), the function G^^ has a pole at cj = ^2{^) 
near which it has the form 

G'\uj ^ e2(q), q, P, k) « ^^^ T, x ^ (12) 



qa 



where 



Map = Sap + JalaP ^ 2FSJa COS —Vp, 



d''z 



jafj 



{2^YV{q,zy 

d'^z cos z" 
{2nrV{q,z)' 

d'^z cos z" cos z^ 



(18) 



w - e2(q) + i5' 



where 



<I>,(p) = (2;q|a 



,„ , a ,„ 10) 

q/2+P q/2-p I / 



{2^Y V{q,z) ' 

and 'D{q, z) is given by Eq. JTol). Then, if determinant of 
M is zero at some w, Eq. (|16[) has a non-trivial solution 
and the vertex has a pole at the corresponding a; and q. 
Otherwise, the residue is zero and there is no bound state 
at the corresponding w and q. 



is a bound state wave function in momentum notation. 



Two-particle Green's function (11) is equal to the vertex 
r(w, q, p, k) multiplied by 



One notes that parameters ( 18 1 depend on D, h and uj 
only in the combination 

h = h-uj/2- D{2S- 1) 



and M{q,D,h,uj) is actually a function of {D,h). Con- 
sidering D and h as independent variables, we note that 
det M{D,h) is a linear function of D. It signifies that 
there exists a unique D value at which det M(I?, ft,) = 
for any h larger than some value ho such that denomi- 
nator (10 1 is positive at h > ho. Then, the two-magnon 
bound state exists ath>ho and the critical field at which 
the two-magnon bound state spectrum becomes gapless 
is higher than the critical field of the one-magnon mode 
(if these critical fields exisiP) . 

The above equations are applicable for consideration 
of two-magnon bound states in strong magnetic field on 
any Bravais lattice with any exchange interactions and 
any S ^ 1. They can be easily extended also to the case 
of anisotropic exchange. We apply the above equations 
now to the particular spin-1 model on a tetragonal lattice 
with different unfrustrated exchange constants along z- 
axis (J^) and in xy plane (J^^) which are parametrized 
as follows: 



Q -3 




FIG. 2. Solid red line (black dots) represents DdO) curve 
(see text). Colored area is the nematic phase stability region. 



J"" = COS 9, 
^y =sin6', 

'TT < 9 ^TT. 



(19) 



The cases oi 9 — and tt correspond to ID AF and FM 
chains, respectively, 9 = ±7r/2 correspond to 2D AFs 
and FMs on the square lattice. Other 9 describe axially 
symmetric 3D magnets. 



General consideration. 



One has Jr, 



2JzC0spz + 2Jxy{cospx + cospy) and we find using the 
above formulas that the bound state at strong field ex- 
ists at any 6* if D < Dc{9). The bound state spec- 
trum is always quadratic near it's minimum p = 0: 
e2(p) « 2{h-hs)+Czpl+C^y{pl+pl). Function Dc (6*) is 
shown in Fig. I2]by red solid line. Behavior of Dc{9) near 
its singularities which correspond to ID and 2D systems 
can be found using (18) with the following result: 



Dc{9 - 0) w -4/3 + const x ^1^1, 

„ /„ , N const 

Dc{9 - ±7r) w const x .J\9^it\. 

The bound state can formally become gapless upon the 
field decreasing if 9 lies in the range — 0.1977r < 9 < 
O.QOStt. In this interval, Bose condensate of bound state 
quasiparticles would appear below the critical field. But 
a spin-flop transition at certain values of D occurs at 
field values which are higher than the critical field of 
the bound state. Detailed consideration of the spin- 
flop transition is out of the aim of the present paper. 
We restrict ourselves by a simple classical result for the 
spin-flop field obtained by minimization of the classical 
energy. According to the classical result, the spin- flop 
transition occurs before condensation of bound states at 
D < Dsfi{9), 9 < and D < Dsf2{9), 9 > tt/2, where 



Dsfi{9) and Dsf2{9) are shown in Fig.|2]by dashed black 
lines. They approach — cxd as 9 —>■ —0 and 9 — >■ tt/2 + 0, 
correspondingly. As a result the area of the nematic 
phase stability turns out to be narrower and it is shown 
in Fig. [2] 

Assuming repulsion interaction between bound states 
quasiparticles near its spectrum minimum (that is rea- 
sonable due to antiferromagnetic interaction) we obtain 
for the nematic phase at h <hg and 9 ^ Q: 



{S-) = {sf) = o 



and 



J2^o{r,-n)S+S+)=2e^ 



l-{S!)^2p^^^-^, 9^n/2 
Inp 



(20) 
(21) 



where $o(^) is the wave function of the bound state ( 14) 
with zero momentum in the spatial representation, p is 
the condensate density, and A is an effective coupling of 
bound state quasiparticles near the spectrum minimum 
(in 3D case). Thus, we obtain a long-range nematic order 
(Eq. (21 1 ) without a conventional magnetic order in xy 
plane (Eq. ( 20 ) ) . Due to the quadratic dispersion of £2 (p) 
we also have the 3D BEG relation^ 



ft,(o)-ft,(T)(xr3/2. 



(22) 



4. ID AF chain (9 = 0). One has in this case 
Jp = 2cosp and all results (saturation field, C^, etc.) 
can be found analytically which however turn out to be 
quite cumbersome. Then, we do not present all of them 
here. As it is mentioned above, bound state exist at 
D < Del — —4/3. The saturation field h^ is shown in 
Fig. Islas a function of D. Bearing in mind the quadratic 
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FIG. 3. Phase diagram of the isolated chain. The Neel and 
Haldane phases are gapped, while the nematic and canted 
antiferromagnetic (CAF) ones are gapless. Inset shows the 
region, where four phases meet. 



dispersion of e2(p) near its minimum and using results 
of ID Bose-gas discussions^ one obtainai^. at /i < /is and 
T = 



1 - (5|) - 2p 



2{hs - h) 



C 



(23) 



(5|+„(i)5|(0)) « (5|)' 



(24) 
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1 



1 

-K \(n + mi)2 (n — iut)^ 

2 /^_,„^^2\ i?2 



Bi 



cos{n{a-aj)ri) 



n^ + u^t 



2+2 



{s^it)) {s-m 



\J\n + iut\ 



(25) 



where p is the density of quasiparticles describing bound 
states, n ^ 1, M = ^ttCzP and B12 are constants. AU 
other correlators are either exponentially small or equal 
to zero exactly. Then, the nematic phase in the isolated 
chain is characterized by no long-range magnetic order 
and quasi-long-range nematic order described by the al- 
gebraic decay of the nematic correlator ( 25 ) . 



To support our analytical results we have performed 
numerical analysis using exact diagonalization applica- 
tion from ALPS software package.^I^ This numerical con- 
sideration is simplified by the fact that the Hamiltonian 
(n]) commutes with the z component of the total spin S^' 
and with the Zeeman term. As a result all the Hamilto- 
nian eigenstates can be classified by eigenvalues M of S^' . 
Let us denote E{M) the minimum energy in each M sec- 
tor dA, h — Q. The ground state energy of a cluster with 
L spins in magnetic field is given by the minimum value 
of E{M) — hM. An important observation is that values 
E{M) — hM at even Mgat — M are smaller than those 



with odd Msat — M, where Mgat — L. This property 
holds for different sizes (L = 6 -;- 16) and both periodic 
(for even L) and open boundary conditions. 

Thus, one can expect that a quasi-condensation takes 
place in the thermodynamical limit of elementary exci- 
tations carrying spin 2. Values of hg obtained numeri- 
cally as a result of analysis of clusters with L = 8 -;- 40 
are in excellent agreement with the corresponding values 
found analytically. Numerical consideration of clu sters 
with L = 8 -h 24 similar to that performed in Refs.'^^IlSl 
confirms also the validity of Eq. (23). For example, one 
obtains for D = -3.5 that 1 - (3^) == a{hs - hf/^ 
where a = 6.18 ± 0.15 and 5 = 2.0041 ± 0.003, that is 
in good agreement with the analytical result 1 — (^S^) w 

6.37-\/iJs — H- We have calculated numerically magne- 
tization near saturation for several D values and found 
results which differ from the theory by a few percents 
only. 

The phase diagram obtained numerically is shown in 
Fig. Is] that is in very good agreement with previous nu- 
merical discussions of the present ID systerrfSHlSl. The 
line of phase transitions from the fully polarized state is 
obtained by consideration of some small-(Msat — M) sec- 
tors of clusters with L = 8 -^ 40. This transition takes 
place to the nematic and CAF phases at Z? < Dd = 
—4/3 and D > Dd, respectively, that is in agreement 
with our theoretical predictions. 

At small h the system goes from the Haldane phase 
to that with Neel long-range order as the easy-axis 
anisotropy rises. The transition point D — Dc2 ~ —0.394 
is found using the extrapolating procedure from L = 
8, 10, 12, 14, 16 as it is explained in section HI.B of Ref.™. 
The value Z?c2 does not depend on h because both Neel 
and Haldane phases have gaps and their ground state 
energy do not depend on h. 

Phase transitions take place upon the field increas- 
ing from Neel and Haldane phases to nematic and CAF 
ones. We observe that the lowest excited states belong 
either to M = 1 or M = 2 sectors in Haldane and Neel 
phases. Then, corresponding critical fields hc{M) can 
be found using the value hc{M, L) = ^ ' ]^ ^ ' ' as 
hc{M,L ^ 00) = hc{M). If /ic(l) < (>)/ic(2) the tran- 
sition takes place to CAF (nematic) phase. For each D 
value, hc{M, L) have been found for L = 8, 10, 12, 14 and 
hc{l), hc{2) have been obtained by an extrapolatiorPSl of 
hc{M, L) to the thermodynamical limit L ^ 00. Wc have 
obtained that hc{l) and hc{2) intersect at Dcs ~ —0.352 
(see Fig.lsl. Two multi-critical points are shown in Fig.ls] 
which are close to each other: 0^2 ~ —0.394 ± 0.01, 
hc2 w 0.817 ± 0.013 and Dc3 » -0.352 ± 0.003, /i^a « 
0.772 ± 0.004. However, the precision of our calculations 
is insufficient to exclude the possibility of these two points 
merging into a single multi-critical one. The boundary 
between the nematic and CAF phases in Fig. [s] (dashed 
line) is taken from Ref.'i^ 

5. Summary and conclusion. We demonstrate in 
the present paper that easy-axis single-ion anisotropy D 
in non- frustrated spin-1 antiferromagnetic systems can 



lead to spin-nematic states in strong magnetic field h. 
General formulas are derived for consideration of transi- 
tions in such systems from the fully polarized phase to 
the nematic one. The region of the nematic phase stabil- 
ity is obtained in the model (fTj) with exchange constants 
(19) parametrized by the value 9 (see Fig. [2]). We discuss 
in detail both analytically and numerically the special 
case oi 9 = that describes AF chain. We obtain the 
phase diagram of the isolated chain on h — D plane (see 



Fig. [3| . We hope that the present study will stimulate 
further theoretical and experimental activity in nematic 
phases discussion both in the considered systems and in 
other ones containing anisotropy. 
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